Calculating Ensemble Averaged Descriptions of Protein Rigidity without Sampling

Previous works have demonstrated that protein rigidity is related to thermodynamic stability, especially under conditions that favor formation of native structure. Mechanical network rigidity properties of a single conformation are efficiently calculated using the integer body-bar Pebble Game (PG) algorithm. However, thermodynamic properties require averaging over many samples from the ensemble of accessible conformations to accurately account for fluctuations in network topology. We have developed a mean field Virtual Pebble Game (VPG) that represents the ensemble of networks by a single effective network. That is, all possible number of distance constraints (or bars) that can form between a pair of rigid bodies is replaced by the average number. The resulting effective network is viewed as having weighted edges, where the weight of an edge quantifies its capacity to absorb degrees of freedom. The VPG is interpreted as a flow problem on this effective network, which eliminates the need to sample. Across a nonredundant dataset of 272 protein structures, we apply the VPG to proteins for the first time. Our results show numerically and visually that the rigidity characterizations of the VPG accurately reflect the ensemble averaged properties. This result positions the VPG as an efficient alternative to understand the mechanical role that chemical interactions play in maintaining protein stability

FlexOracle: predicting flexible hinges by identification of stable domains

<p>Abstract</p> <p>Background</p> <p>Protein motions play an essential role in catalysis and protein-ligand interactions, but are difficult to observe directly. A substantial fraction of protein motions involve hinge bending. For these proteins, the accurate identification of flexible hinges connecting rigid domains would provide significant insight into motion. Programs such as GNM and FIRST have made global flexibility predictions available at low computational cost, but are not designed specifically for finding hinge points.</p> <p>Results</p> <p>Here we present the novel FlexOracle hinge prediction approach based on the ideas that energetic interactions are stronger <it>within </it>structural domains than <it>between </it>them, and that fragments generated by cleaving the protein at the hinge site are independently stable. We implement this as a tool within the Database of Macromolecular Motions, MolMovDB.org. For a given structure, we generate pairs of fragments based on scanning all possible cleavage points on the protein chain, compute the energy of the fragments compared with the undivided protein, and predict hinges where this quantity is minimal. We present three specific implementations of this approach. In the first, we consider only pairs of fragments generated by cutting at a <it>single </it>location on the protein chain and then use a standard molecular mechanics force field to calculate the enthalpies of the two fragments. In the second, we generate fragments in the same way but instead compute their free energies using a knowledge based force field. In the third, we generate fragment pairs by cutting at <it>two </it>points on the protein chain and then calculate their free energies.</p> <p>Conclusion</p> <p>Quantitative results demonstrate our method's ability to predict known hinges from the Database of Macromolecular Motions.</p

Changes in Lysozyme Flexibility upon Mutation Are Frequent, Large and Long-Ranged

We investigate changes in human c-type lysozyme flexibility upon mutation via a Distance Constraint Model, which gives a statistical mechanical treatment of network rigidity. Specifically, two dynamical metrics are tracked. Changes in flexibility index quantify differences within backbone flexibility, whereas changes in the cooperativity correlation quantify differences within pairwise mechanical couplings. Regardless of metric, the same general conclusions are drawn. That is, small structural perturbations introduced by single point mutations have a frequent and pronounced affect on lysozyme flexibility that can extend over long distances. Specifically, an appreciable change occurs in backbone flexibility for 48% of the residues, and a change in cooperativity occurs in 42% of residue pairs. The average distance from mutation to a site with a change in flexibility is 17–20 Å. Interestingly, the frequency and scale of the changes within single point mutant structures are generally larger than those observed in the hen egg white lysozyme (HEWL) ortholog, which shares 61% sequence identity with human lysozyme. For example, point mutations often lead to substantial flexibility increases within the β-subdomain, which is consistent with experimental results indicating that it is the nucleation site for amyloid formation. However, β-subdomain flexibility within the human and HEWL orthologs is more similar despite the lowered sequence identity. These results suggest compensating mutations in HEWL reestablish desired properties

Common Variants at 10 Genomic Loci Influence Hemoglobin A(1C) Levels via Glycemic and Nonglycemic Pathways

OBJECTIVE Glycated hemoglobin (HbA1c), used to monitor and diagnose diabetes, is influenced by average glycemia over a 2- to 3-month period. Genetic factors affecting expression, turnover, and abnormal glycation of hemoglobin could also be associated with increased levels of HbA1c. We aimed to identify such genetic factors and investigate the extent to which they influence diabetes classification based on HbA1c levels. RESEARCH DESIGN AND METHODS We studied associations with HbA1c in up to 46,368 nondiabetic adults of European descent from 23 genome-wide association studies (GWAS) and 8 cohorts with de novo genotyped single nucleotide polymorphisms (SNPs). We combined studies using inverse-variance meta-analysis and tested mediation by glycemia using conditional analyses. We estimated the global effect of HbA1c loci using a multilocus risk score, and used net reclassification to estimate genetic effects on diabetes screening. RESULTS Ten loci reached genome-wide significant association with HbA1c, including six new loci near FN3K (lead SNP/P value, rs1046896/P = 1.6 × 10−26), HFE (rs1800562/P = 2.6 × 10−20), TMPRSS6 (rs855791/P = 2.7 × 10−14), ANK1 (rs4737009/P = 6.1 × 10−12), SPTA1 (rs2779116/P = 2.8 × 10−9) and ATP11A/TUBGCP3 (rs7998202/P = 5.2 × 10−9), and four known HbA1c loci: HK1 (rs16926246/P = 3.1 × 10−54), MTNR1B (rs1387153/P = 4.0 × 10−11), GCK (rs1799884/P = 1.5 × 10−20) and G6PC2/ABCB11 (rs552976/P = 8.2 × 10−18). We show that associations with HbA1c are partly a function of hyperglycemia associated with 3 of the 10 loci (GCK, G6PC2 and MTNR1B). The seven nonglycemic loci accounted for a 0.19 (% HbA1c) difference between the extreme 10% tails of the risk score, and would reclassify ∼2% of a general white population screened for diabetes with HbA1c. CONCLUSIONS GWAS identified 10 genetic loci reproducibly associated with HbA1c. Six are novel and seven map to loci where rarer variants cause hereditary anemias and iron storage disorders. Common variants at these loci likely influence HbA1c levels via erythrocyte biology, and confer a small but detectable reclassification of diabetes diagnosis by HbA1c

How well do blood folate concentrations predict dietary folate intakes in a sample of Canadian lactating women exposed to high levels of folate? An observational study

<p>Abstract</p> <p>Background</p> <p>In 1998, mandatory folic acid fortification of white flour and select cereal grain products was implemented in Canada with the intention to increase dietary folate intakes of reproducing women. Folic acid fortification has produced a dramatic increase in blood folate concentrations among reproductive age women, and a reduction in neural tube defect (NTD)-affected pregnancies. In response to improved blood folate concentrations, many health care professionals are asking whether a folic acid supplement is necessary for NTD prevention among women with high blood folate values, and how reliably high RBC folate concentrations predict folate intakes shown in randomized controlled trials to be protective against NTDs. The objective of this study was to determine how predictive blood folate concentrations and folate intakes are of each other in a sample of well-educated lactating Canadian women exposed to high levels of synthetic folate.</p> <p>Methods</p> <p>The relationship between blood folate concentrations and dietary folate intakes, determined by weighed food records, were assessed in a sample of predominantly university-educated lactating women (32 ± 4 yr) at 4-(n = 53) and 16-wk postpartum (n = 55).</p> <p>Results</p> <p>Median blood folate concentrations of all participants were well above plasma and RBC folate cut-off levels indicative of deficiency (6.7 and 317 nmol/L, respectively) and all, except for 2 subjects, were above the cut-off for NTD-risk reduction (>906 nmol/L). Only modest associations existed between total folate intakes and plasma (r = 0.46, <it>P </it>< 0.001) and RBC (r = 0.36, <it>P </it>< 0.01) folate concentrations at 16-wk postpartum. Plasma and RBC folate values at 16-wk postpartum correctly identified the quartile of folate intake of only 26 of 55 (47%) and 18 of 55 (33%) of subjects, respectively. The mean RBC folate concentration of women consuming 151–410 μg/d of synthetic folate (2^ndquartile of intake) did not differ from that of women consuming >410 μg/d (3^rdand 4^thquartile).</p> <p>Conclusion</p> <p>Folate intakes, estimated by food composition tables, and blood folate concentrations are not predictive of each other in Canadian lactating women exposed to high levels of folate. Synthetic intakes > 151–410 μg/d in these women produced little additional benefit in terms of maximizing RBC content. More studies are needed to examine the relationship between blood folate concentration and NTD risk. Until data from such studies are available, women planning a pregnancy should continue to consume a daily folic acid supplement of 400 μg.</p

Free backbone carbonyls mediate rhodopsin activation

Conserved prolines in the transmembrane helices of G-protein-coupled receptors (GPCRs) are often considered to function as hinges that divide the helix into two segments capable of independent motion. Depending on their potential to hydrogen-bond, the free C=O groups associated with these prolines can facilitate conformational flexibility, conformational switching or stabilization of the receptor structure. To address the role of conserved prolines in family A GPCRs through solid-state NMR spectroscopy, we focus on bovine rhodopsin, a GPCR in the visual receptor subfamily. The free backbone C=O groups on helices H5 and H7 stabilize the inactive rhodopsin structure through hydrogen-bonds to residues on adjacent helices. In response to light-induced isomerization of the retinal chromophore, hydrogen-bonding interactions involving these C=O groups are released, thus facilitating repacking of H5 and H7 onto the transmembrane core of the receptor. These results provide insights into the multiple structural and functional roles of prolines in membrane proteins

Probing the Flexibility of Large Conformational Changes in Protein Structures through Local Perturbations

Protein conformational changes and dynamic behavior are fundamental for such processes as catalysis, regulation, and substrate recognition. Although protein dynamics have been successfully explored in computer simulation, there is an intermediate-scale of motions that has proven difficult to simulate—the motion of individual segments or domains that move independently of the body the protein. Here, we introduce a molecular-dynamics perturbation method, the Rotamerically Induced Perturbation (RIP), which can generate large, coherent motions of structural elements in picoseconds by applying large torsional perturbations to individual sidechains. Despite the large-scale motions, secondary structure elements remain intact without the need for applying backbone positional restraints. Owing to its computational efficiency, RIP can be applied to every residue in a protein, producing a global map of deformability. This map is remarkably sparse, with the dominant sites of deformation generally found on the protein surface. The global map can be used to identify loops and helices that are less tightly bound to the protein and thus are likely sites of dynamic modulation that may have important functional consequences. Additionally, they identify individual residues that have the potential to drive large-scale coherent conformational change. Applying RIP to two well-studied proteins, Dihdydrofolate Reductase and Triosephosphate Isomerase, which possess functionally-relevant mobile loops that fluctuate on the microsecond/millisecond timescale, the RIP deformation map identifies and recapitulates the flexibility of these elements. In contrast, the RIP deformation map of α-lytic protease, a kinetically stable protein, results in a map with no significant deformations. In the N-terminal domain of HSP90, the RIP deformation map clearly identifies the ligand-binding lid as a highly flexible region capable of large conformational changes. In the Estrogen Receptor ligand-binding domain, the RIP deformation map is quite sparse except for one large conformational change involving Helix-12, which is the structural element that allosterically links ligand binding to receptor activation. RIP analysis has the potential to discover sites of functional conformational changes and the linchpin residues critical in determining these conformational states